Degeneracy in Density Functional Theory: Topology in v- and n-Space 
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This paper clarifies the topology of the mapping between v- and n-space in fermionic systems. 
Density manifolds corresponding to degeneracies g — 1 and g > 1 are shown to have the same 
mathematical measure: every density near a g-ensemble-w-representable (<?-VR) n(r) is also fir-VR 
(except "boundary densities" of lower measure). The role of symmetry and the connection between 
T = and T = + are discussed. A lattice model and the Be-series are used as illustrations. 
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Density functional theory (DFT) of electrons in an ex- 
ternal potential v(r) uses the ground-state density n(r) 
as basic variable [jl], ^]. n{r) uniquely determines u(r) 
(apart from a constant); conversely, for a given number, 
N, of electrons, v(r) obviously determines n(r) uniquely 
if and only if the ground state is non-degenerate. These 
facts point to the interesting role of degeneracy in the 
mapping i;(r)<S=> n(r). 

In 1982-3, Levy || and Lieb Q discovered exam- 
ples of well-behaved density functions, n(r), which could 
not be reproduced as non-degenerate, non-interacting 
ground-state densities of any v(r), and were called non- 
w-representable (non-VR). In these examples, the n(r) 
could be reproduced as weighted averages of densities of 
degenerate ground states corresponding to a v(r). 

Soon afterwards, Kohn || adopted a lattice version 
of the Schrodinger equation, r — ► n (I = 1, . . . , M), in 
which v(r) and n(r) become vi and ni and can be viewed 
as M'-dimensional (M'-D) vectors: v = (v\, . . . , Vm'), 
with M' = M — 1, ujf = 0; and n = (ni, . . . , ?Hf/), with 
riM = N — [nx + . . . + riM') @- The M lattice points 
are enclosed by boundary points on which v = +oo and, 
accordingly, all wave functions vanish. It was shown in 
H that in the M'-D n-space there are finite M'-D regions 
in which all "points" (densities) are VR. 

In 1985, Chayes et al. (CCR) (tJ, using a similar lattice 
model, proved the following important result: Any well- 
behaved n (0 < ni < 1) can be uniquely represented as 
weighted average of g degenerate ground-state densities 
associated with a v. We call the special case g = 1 1- 
VR, and the general case (g > 1) g-VR. Thus, n-space is 
filled by manifolds S 9 , in each of which n is g-VR. Each 
S 9 maps on a corresponding manifold Q 9 in ti-space. 

This paper aims at an understanding of the topolo- 
gies of the regions S 9 and Q 9 and of the nature of the 
boundaries between them. This has mathematical and 
physical significance, e.g. in the search for self-consistent 
solutions of Kohn-Sham (KS) equations M. We shall see 
that g-VR densities with g > 1 are not mathematically 
exceptional and, in particular, do not depend on symme- 
tries of the potential. As illustrations we shall present a 
finite lattice model and the Be-series for Z — > oo M. We 
also discuss the relationship between T = ensembles 



and T = + thermal ensembles. We limit ourselves to 
systems that are everywhere nonmagnetic. 

Interior regions in n-space. We first recapitulate 
Ref. H for 1-VR densities, and then generalize to g-VR. 

a) g = 1. Since there is a non-degenerate non- 
interacting ground state with a finite gap, corresponding 
to some v, we can use non-degenerate perturbation the- 
ory to calculate the first-order density change, n;, due 
to a weak perturbing potential 1 vf. 
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Since, according to yj], the density n of a non-degenerate 
ground state uniquely determines v (no arbitrary con- 
stant since vm = 0), the homogeneous set of equations 
corresponding to ([!]) has no solution other than x vv = 0, 
so that (111) can be inverted: 
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Thus, any first-order change of n preserving N produces 
a density which is also 1-VR. This is the case in the entire 
M'-D neighborhood surrounding the point n, in which A 
remains positive. 

b) g > 1 0. We denote by . . . , * g a set of g 
orthogonal degenerate ground-state wave functions cor- 
responding to v 9 , with common energy E and finite gap 
A. The remaining eigenfunctions are ^> g+ i,.... What 
conditions must be imposed on infinitesimal [O(A)] po- 
tential changes v which will preserve this degeneracy? 
In the space of the \t j , the perturbed Hamiltonian matrix 
is 



E i 5 il + 1 V l 
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Ei=E (i= 1,.. .,g), Ei > E + A (i > g + 1), and x V l} 



JZftiii^vilj). Because of the finite A, the off-diagonal 
Vij with i < g and j > g + 1, or vice versa, can be 
removed to all orders in A by orthogonal transformations 
leading to the decoupled g x g block Hamiltonian 



l Vtj + 0(X 2 ) 
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Diagonalization by an orthogonal transformation = 
tij+0(X 2 ) and the requirement that the eigenvalues 
remain degenerate gives a transformed Hamiltonian 



fly = (T-^T) = (E + 1 e)S ij + 0(X 2 ) 



(5) 



On inverting the transformation and noting that every 
matrix commutes with Sjj, one finds Hij = Hij + 0(X 2 ). 
Comparison with Eq. (Q) leads to \ (g — 1) (g + 2) con- 
ditions of equal diagonal and vanishing off-diagonal Vy. 
Thus, every point v 9 with a 5-fold degenerate ground 
state and finite A is embedded in a manifold in u-space 
of dimension D 9 = \M' — \{g — l)(g + 2)] , in which the 
degeneracy g and a finite A are preserved. Ground-state 
degeneracies in w-space are "rare" in the above sense. 
Each v 9 gives rise to a set of ensemble densities, 
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< w j < 1, 



a 



= 1 



(6) 



with rij — rij (R 9 ^j), where R 9 is a g-D orthogonal 
transformation [^g(g — 1) parameters] and the Wj are 
normalized weights {g — 1 parameters). Thus, all ft' in 
a finite A/'-D neighborhood enclosing a g-VR ft, defined 
by A(n') > and Wj > 0, are also <?-VR: Degeneracy in 
n-space is not "rare" . 

Boundary surfaces in n-space. Except for bound- 
aries of D < M' — 1, n-space is completely filled by M'-D 
regions S 9 , with degenerate ground-state levels and a fi- 
nite positive gap to the nearest excited state. 

Apart from "corners" (D < M' — 2), interior regions 
S 9 in n-space are bounded by (M' — 1)-D internal and 
external "surfaces" S s , see Fig. [|. Each internal E 9 
separates two interior regions S 9 and S 9+1 . A point ft 
on such a S 9 corresponds to a (g + l)-fold degenerate 
ground state but with one of the g + 1 weights equal to 
zero. On the g-side of Y, 9 there is a <?-fold degeneracy and 
a small gap to the (g + l)st state opens up from 0. On 
the (g + l)-side there is a (<? + l)-fold degeneracy and one 
of the g + 1 weights starts from and becomes positive. 

External boundaries correspond to either one of the 
n\ becoming zero, reflecting particle conservation; or one 
of the ni becoming 1, reflecting the Pauli principle. For 
examples of both kinds, see the lattice model below. 

Imposed symmetry conditions. Consider a Hamil- 
tonian, invariant under a group G of the lattice with ir- 
reducible representations h with D = g^. Let the ground 
state have total degeneracy g — J2h m hgh due to rrih oc- 
currences of h. If a small perturbation H ' respecting the 
group G is imposed, how many further conditions must 
be met by H' so that the degeneracy is maintained? 

We have shown above that the perturbation must have 
vanishing off-diagonal matrix elements (MEs) between all 
degenerate ground states, and equal diagonal MEs. 

Off-diagonal MEs. The symmetry of H' assures imme- 
diately that all off-diagonal MEs between h ^ h 1 vanish 
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FIG. 1: Two interior regions in the A/'-D n-space, S 9 and 
S 9+1 , separated by a (M' - 1)-D boundary E 9 . 



JTT| . Off-diagonal MEs between two states of the same h 
involve, by the Wigner-Eckart Theorem, only one inde- 
pendent parameter. The total number of independent pa- 
rameters characterizing all off-diagonal MEs between the 
degenerate ground states is therefore '^2 h mh(mf l — l)/2. 

Diagonal MEs. The requirement that all diagonal MEs 
of H' are equal is automatically satisfied for states be- 
longing to the same h, and involves ^2 h mh — 1 conditions 
for states belonging to different h. 

Thus, to ensure that the g-fold degenerate ground state 
maintains its full degeneracy requires the imposition of 
Y^hl m h( m h — l)/2 + mh] ~ 1 conditions. 

Next we ask: for a ground state with g = X)/i m '»5'»' 
how many free parameters are there to create ensemble 
densities invariant under group G of the Hamiltonian? 

The most general ensemble density obtainable from g 
degenerate ground states ^1, . . . , ^ g is YjL=i w h ^O^*)) 
= R 9 '$>[_ 1 . However, the must be chosen in a re- 
stricted way to assure that the density has the full sym- 
metry of the Hamiltonian: Let l ^.' Uh (k — 1, . . . , gu) be a 
set of orthogonal eigenfunctions, constituting a basis for 
a given h of the group G, with occurrence 1/^ = 1,..., m^, 
such that all ty?^'"' 1 transform in the same way under G. 



The most general density of the required symmetry is 
n = Y J wh,Uh H ( RUh( Pk' /h ) ■ ( 7 ) 



The total number of free parameters is given by those in 
the rotation R Uh plus the number of independent weights 
w h '" h , i.e. Y^hl mh ( rnh ~ + — 1. We observe that 
the total number of constraints on H' is equal to the total 
number of free parameters in ft. Therefore the dimension- 
ality of the density invariant under G associated with any 
set of degenerate states is the same, independent of G. 

Lattice Example. We illustrate our topological re- 
sults in an example including symmetry requirements in 
v- and n-space. Consider the singlet ground state of 4 
non-interacting electrons, or, equivalently, 2 spinless non- 
interacting fermions, whose wave functions (pj are defined 
on a 13-point fee lattice (see Fig. ||), and vanish on all 
exterior points. The discrete Schrodinger equation is 
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FIG. 2: 13-point fee lattice with dh symmetry: Vi = Vi 
(V 2 ) on the sites I = 1-8 (9-12), and vi 3 = 0. 



where, for fee lattices, V 2 /; = [^ fc /& — 12/)]/a 2 , with fc 
running over the nearest neighbors of lattice site I. Ener- 
gies and potentials will be measured in units of 2ma 2 /h 2 . 

We consider cases of C^h symmetry: vi = V\, I = 1-8; 
Vi = V%, I = 9-12; V13 = 0. The lattice wave functions 



belong to the following representations of C±h |ll|] : A g (s- 
like), E u and A u (p-like), E g and B g (d-like), and B u (/- 
like). A g occurs three times (lA g , 2A g , 3A g ), E u occurs 
twice (1E U , 2E U ), all other representations occur only 
once. The eigenfunctions and energies required solutions 
of 3-, 2- and 1-D secular equations. 

Consider the ground state of 2 non-interacting spinless 
fcrmions. The lowest level is always lA g (ls-like). The 
next level depends on V\ and V2, see Fig. ||. The (Vl, V2) 
plane divides into three distinct, infinitely extended re- 
gions (shown in yellow), in each of which the eigenfunc- 
tion of the second level belongs to a particular irreducible 
representation of C4/,: Region la, 1A U (2p z -like); region 
lb, 2A g (2s-like); and region lc, 1E U (2p X7 2p y -\ike). In- 
side each region, the ground state is non-degenerate (not 
counting degeneracies dictated by symmetry, such as in 
the 2-D representation 1E U ). 

On the boundary lines 2a, 2b and 2c (shown in blue), 
crossings of levels belonging to different representations 
of C±h occur, leading to two-fold (accidental) degen- 
eracies. Along 2a, the potential has cubic symmetry 
(Vi = V2), and the two lowest p-like levels (1A U , 1E U ) 
coincide. On 2b, the 2A g and 1A U levels cross, and on 
2c, the 2A g and 1E U levels cross. All three boundary 
lines meet in a single point, "3" , at V\ = V2 = 8. At this 
point only, the ground state has a three-fold accidental 
degeneracy. Note that the degeneracies on 2b and 2c are 
not due to additional potential symmetries. 

As shown above, a ground state with restricted symme- 
try maintains its degeneracy if Ylhl m h( m h~ 1)/2+to/J— 1 
conditions are met. In the present example, = 1 al- 
ways. Regions 2a, 2b, 2c are therefore lines (1 condition), 
region 3 is a point (2 conditions) in the (Vi, V2) -plane. 

We now discuss all ft with C/^h symmetry. Like v, they 
are fully characterized by two values, N% and A2, the 
total density at points 1-8 and 9-12. 7113 = 2 — Ni — N2- 
N% and are confined within a stripe in the (Nx,N2) 
plane, with an "inner" boundary, with n.13 = 1, reflecting 
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FIG. 3: Top: ground states of 2 non-interacting spinless 
fermions on the lattice of Fig. ^, with potentials of dh sym- 
metry. Bottom: associated region of allowed densities. 



the Pauli principle (Ni + N 2 < 1 would imply 77.13 > 
1), and an "outer" boundary, with riia = 0, reflecting 
particle conservation (Ni + N 2 < 2). 

The manifold of allowed densities in the (N\, N2) plane 
is subdivided into different regions, corresponding to 
those of the (Vi, VO-plane. All of the regions in the 
(JVijiVa) plane have D = 2. This is in line with our 
general topological results for g-VK densities. 

Fig. ^ illustrates the highest occupied state along 
paths in the (VijVz) and (A^iVa) planes (see dotted 
lines in Fig. ||top and bottom): 1A U with weight 1 (re- 
gion la), 1E U with weight 1 (equal- weight combination of 
2p Xl 2p y -like states; region lc), or a linear combination of 
the degenerate 1A U and 1E U (region 2a), with fractional 
weights. 

Physical Example: the Be-series. Be has a nomi- 
nal electronic configuration (ls) 2 (2s) 2 and a term value 
1 S. Without interactions, 2s and 2p are degenerate. 
The Be-series consists of neutral Be (Z = 4)|l2|, and 
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FIG. 4: Energies and ensemble weights of the next two levels 
above lA g , along paths in the (Vi,V2) and (Ni,N2) planes 
(dotted lines in Fig. 0). 



the 4-electron ions with Z = 5, 6, ... , oo. By re-scaling 
(f = Zr/4), this is equivalent to a series of Be-like atoms 
with Z = 4 but with weakened interaction strength 
e 2 = e 2 (A/Z). For Z — > oo, only the nearly degen- 
erate configurations (ls) 2 (2s) 2 and (ls) 2 (2p) 2 , (2p) 2 — 
3-1/2 J2^ = i(2pj) 2 , need to be considered. This leads to 
a 2 x 2 secular equation all of whose MEs are exactly 
known mM . The limiting density is explicitly 4-VR: 



n (r) 



Z\2 



<^i s (f) 2 + 0.949^ 2s (f) 



0.051^ s 



(9) 

where the y>'s are scaled hydrogenic wave functions (for 
Z = 4). The configurations (ls) 2 (2s) 2 and (ls) 2 (2p) 2 
have amplitudes 0.974 and 0.225. (Unlike in CI expan- 
sions jbj , KS densities consist strictly of a finite number 
of orbital densities for all Z) . 

From our previous discussion of boundary surfaces, we 
observe that the density (^|) is not on a boundary, since 
there is no unoccupied degenerate state. Under these 
conditions, any density n(r) sufficiently near no(r), in 
general without any symmetry, is also 4-VR. In particu- 
lar, there must be a finite range in 1/Z of the Be-series 
where densities are 4-VR. This is borne out by numerical 
calculations 

Density ensembles at T = and + . Let v 9 be 

a potential with g degenerate non-interacting, gapped 
ground states, $j, generating the i(g — l)(g + 2)-D set 
of ensemble densities (0). Now introduce a weak per- 



turbing potential V of order A which, in general, splits 
the degeneracy and leads to eigenstates and eigenvalues 



0(A); 



E 3 = E 



x Ej, where % = 0(A), 



and the corresponding densities are fij. 

For a given low temperature, (3~ l <C A, the canonical 
density corresponding to v + v is given by 



= Z 



9 

_1 E e 
3=1 



^ n{R s ^j) + 0(A) + (e- fiA ) , (10) 



where Z = Y^j e ^ ^ . This density is equal, to lead- 
ing order, to the density ft of Eq. (||), if 1 u is chosen so 
that R 9 <S>j = *j and Z^e'^^i = Wj. As 0A -> 00, 
1 w ~ — > and n' — » n. This establishes the cor- 
respondence between ensemble densities at T = and 
canonical densities at low temperatures, (3 A 3> 1. 

Concluding remarks. The problem of w-represen- 
tability || ^] had cast a shadow of uncertainty over DFT, 
which was very largely clarified in Q|. Here, building on 
CCR, we clarify the topology of the i?(r)-n(r) mapping. 
The Be-series and a lattice model are used for illustra- 
tion. The transition from a dense discrete lattice, rj, to 
a continuous variable, r, does not appear to offer difficul- 
ties. An open issue is extension to spin magnetism. 
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